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We study the effect on the cosmic microwave background (CMB) anisotropy and iarge scafe 
structure (LSS) power spectrum of a scattering interaction between cold dark matter and baryons. 
This scattering alters the CMB anisotropy and LSS spectrum through momentum transfer between 
the cold dark matter particles and the baryons. We find that current CMB observations can put 
an upper limit on the scattering cross section which is comparable with or slightly stronger than 
previous disk heating constraints at masses greater than 1 GeV, and much stronger at smaller 
masses. When large-scale structure constraints are added to the CMB limits, our constraint is more 
stringent than this previous limit at all masses. In particular, a dark matter-baryon scattering 
cross section comparable to the "Spergel-Steinhardt" cross section is ruled out for dark matter mass 
greater than 1 GeV. 

PACS numbers: 98.80.-k,12.60.-i,14.80.-j 



I. INTRODUCTION 

It is now generally accepted that the energy density 
of the universe includes a substantial fraction of cold 
dark matter (CDM), constituting roughly 0.3 — 0.4 of the 
closure density. It is normally assumed that the CDM 
particle does not interact with ordinary baryonic matter 
or with photons; hence the "dark" in cold dark matter. 
However, a small scattering cross section cannot be en- 
tirely ruled out, only constrained by observations. 

Several different models are possible for interactions 
between ordinary matter and CDM. One possibility is 
strongly interacting dark matter (also known as strongly 
interacting massive particles, or SIMPs), in which the 
CDM particle couples to baryons, but not to photons 
or electrons. Recently, to solve the small scale problem 
in CDM models, strongly self-interacting particles were 
suggested [Q, |), with 



ok a/cm 
8 x 10" 25 < —!— — < 10" 23 . 



m/GeV 



(1) 



is narrower 



The range quoted here (from Ref. is 
than originally proposed in Q| because of additional con- 
straints. While the original Spergel-Steinhardt cross sec- 
tion applies only to dark matter self-interaction, if such 
an interaction is mediated by the strong force, then the 
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dark matter might also interact strongly with baryons M. 
Examples of such dark matter includes, e.g., quark-gluino 
bound states Jp) , strangelets Q , gauge singlet mesons [|| , 
and Q-balls jfjj. Limits on such interactions were inves- 
tigated by Starkman ct al. 0, and recently reexamined 
in references P, |q, pj fiof. A second possibility is to cou- 
ple the dark matter particle electromagnetically, e.g., by 
giving it a tiny electromagnetic charge (see, e.g., refer- 
ence |nj and references therein); such a CDM particle 
will scatter primarily off of photons and free electrons. 

Any interactions of this sort will affect the spectrum 
of the cosmic microwave background (CMB) fluctuations, 
since they act to transfer momentum from the dark mat- 
ter to the baryon-photon fluid, although for WIMPs this 
effect is negligible Jl2], Current CMB observations 
can place constraints on such models. The case of di- 
rect interactions between the CDM particle and photons 
has already been examined pH , so we do not consider it 
here. While the model discussed in reference |l4| is not 
identical to the case of electromagnetically charged dark 
matter, since it includes photon-CDM interactions but 
not electron-CDM interactions, we expect it to be quali- 
tatively similar to the charged dark matter model. There- 
fore, we consider only the case of strongly-interacting 
dark matter, in which the dark matter particle couples 
to baryons (only). 

In the next section, we give the perturbation equations 
for CDM which couples to baryons, and we show the mat- 
ter transfer function and CMB power spectra for some 
representative cases. In Sec. 3, we compare our calcula- 
tions with the CMB observations to derive constraints on 
the cross section for dark matter-baryon scattering, and 
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we compare to previously-derived limits. Our conclusions 
are summarized in Sect. 4. 



II. FLUCTUATION GROWTH WITH 
STRONGLY-INTERACTING DARK MATTER 

A. The Basic Modified Equations 

We begin by assuming a scattering cross section a be- 
tween the CDM particle of mass m c and protons with 
mass m&, and we derive the changes to the perturbation 
equations produced by this scattering interaction. Our 
starting point is the set of perturbation equations in ref- 
erence |5), which form the basis for CMBFAST @. We 
will not reproduce here all of the perturbation equations, 
but simply give those equations which are modified. 

Following reference Jlq] , we work in the synchronous 
gauge. However, because the CDM is interacting it does 
not provide a natural way of defining the synchronous 
coordinates, as opposed to standard CDM |l5| . Let S c 
and 5b be the density fluctuation (= Sp/p) for the cold 
dark matter and the baryons, respectively, and let 9 C and 
9b be the corresponding velocity divergences. Then for 
the cold dark matter, we have 
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For baryons the corresponding terms are 
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In these equations c c and Cb are the adiabiatic sound 
speed of the CDM and baryons, respectively, w rc i is 
the average relative velocity between the CDM and the 
baryons, and R c ,Rb,R-y are the momentum transfer 
terms. The sound speeds are given by 



(k B T b /p) 1 



{k B T c /m c 
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the average relative velocity by 



1 3(m c T b + m b T c 
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In calculating average velocities we have assumed ther- 
mal Maxwcll-Boltzmann distributions for both CDM and 
baryons. The temperature evolution of the dark matter 
and baryons is given by 



Tb + 2-n = -K 7 (T b — T 7 ) — K c (T b — T c ), (12) 
a 



T c + 2-T c = -K b (T c -T b ). 
a 



(13) 



It is possible to obtain analytical solutions for the 
temperature evolution equations, but they are compli- 
cated and cumbersome to use. Instead, we use a semi- 
implicit scheme to integrate these differential equations. 
We tested this method in a few cases, and found that it 
produces solutions which agree with the analytical solu- 
tions within 0.2%, and the effect on the CMB spectrum 
is negligible. 



B. The Effect of Primordial Helium 

This set of equations ignores the fact that the baryonic 
fluid consists of two species: hydrogen and helium-4, in a 
roughly 3 to 1 ratio by mass. To account for the presence 
of helium as well as hydrogen, we make the following 
replacements: 



mbnb<TV rc \ _^ ninTLn <th vh + m Hc nHcO'HcfHc 

nib + Tn c win + m c mne + w c ' 

— > (1 — Y) H Y, (15) 

nib + m c mn + m c toho + m c 



where Y is the primordial mass fraction of helium-4 (we locities for the hydrogen and helium relative to the dark 
take Y = 0.24 throughout), and we define separate ve- matter: 
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In our calculations, <r, the scattering cross section be- 
tween the CDM particle and the proton, is a free pa- 
rameter which we will constrain from the observations. 



However, we need to make some sort of assumption re- 
garding a He, the scattering cross section between CDM 
and helium nuclei. For coherent scattering, we expect 
<THe = 16cr, while a spin-dependent cross section would 
yield an e = 0. For now, we will examine the intermediate 
case of incoherent scattering, for which ane — 4<r. Later, 
in Sec. 3, wc will extend our calculation to consider both 
coherent scattering and spin-dependent scattering. 

Taking mn c ~ 4mn, and n b = ph/mn, the above sub- 
stitution can be achieved by simply replacing everywhere 
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C. Tight-Coupling Approximation photons. In the tight-coupling regime, the temperature 

evolution is given by 

Following reference jl5) , we use a separate set of equa- 
tions in the limit of tight coupling between baryons and 



f b = —2—T b + —-^——an e <TT(Ty — T b ) + K C (T C — T b ) (20) 
a 3 m e p b 

f c = -2-T c + K b (T b -T c ) (21) 
a 



I 

where /i is the mean molecular weight. During tight cou- Setting R = 4/J 7 /3ph, and r 7 = (an e <Tr) _1 , we have 
pling, T c ~ T b ~ T 7 , f b ~ -fT b , and u rc i = -sfwrel- 
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D. Applicable range of our calculation 



with density parameter Q m . We have also assumed that 



In the calculation described above, we have assumed 
that the dark matter is made of a single type of particle, 
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it is a "cold dark matter" particle, i.e. during the whole 
range of our calculation it is non-relativistic, which trans- 
lates to m > 1 MeV. We have also assumed that during 
each scattering only one baryon interacts with the dark 
matter particle. If the interaction length of the particle 
is a (a = Ana 2 ), then we require that this interaction 
length be less than the inter-particle spacing: 



a < n 



-1/3 



(26) 



where rib ~ 10" 5 (1 + z) 3 ft,Qf2bCm~ 3 . Processes which af- 
fect CMB anisotropy happen at redshift less than 10 5 , 
so we can put a limit on the applicable cross section, 
neglecting factors of order unity: 



10" 6 cm 2 . 



(27) 



Our calculation must be modified for a > er max to take 
into account multi-particle scattering. 



E. The matter and CMB power spectra 

We have modified the CMBFAST code |l§ to include 
the effects discussed in the previous sections. For initial 
conditions, we assume that the dark matter and baryons 
started tightly coupled with the same temperature. Ten- 
sor modes have been ignored, and we consider only flat 
geometries (f2 = 1). 

In Fig. we show the CMB fluctuation spectra and 
the corresponding square of the matter transfer function 
for some representative cases. The coupling of the dark 
matter to the baryons damps the matter power spectrum 
on small scales, as can be seen in Fig. 0(b). 

One might imagine that the effect on the CMB of cou- 
pling the baryons to the dark matter would be equivalent 
to a standard model with a larger value of f^. This is 
not the case, as can be seen in Fig. 0(a). Although both 
models produce an increase in the amplitude of the first 
acoustic peak, the large f2f, model produces an increase 
in the amplitude of the third peak, while the strongly- 
interacting dark matter model yields a decrease in the 
third peak amplitude (and in all of the other peaks as 
well). 

One can understand the difference between the CMB 
spectra in the following way: The baryon-photon oscilla- 
tions are damped oscillations, with a damping rate (vis- 
cosity) inversely proportional to the baryon-photon in- 
teraction rate. Interacting CDM with intermediate cross 
section amounts to adding "baryons" but simultaneously 
also increasing the viscosity of the plasma. Therefore 
the acoustic oscillations damp quickly and all peaks but 
the first are suppressed. Scales around the first acous- 
tic peak have only undergone part of an oscillation and 
therefore have not been damped significantly. On such 
large scales the increase in oscillation amplitude due to 
added "baryons" is more important and the overall power 
spectrum increases. 
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FIG. 1: (a) CMB power spectra for a few sample cases. The 
solid (black) curve is a fiducial ACDM model with parameters 
£7 = 1, fl b h 2 = 0.018, Q m h 2 = 0.18, H = 75 kms" 1 Mpc" 1 , 
n — 1, and without dark matter-baryon interaction. The dot- 
ted (blue) curve is for a dark matter mass of 1 GeV, with 
a dark matter-baryon cross section a = 3 x 10 -24 cm 2 . The 
dashed (red) curve is for a model without dark matter-baryon 
interaction, but with increased baryon density, fi;, = 0.05. 
All spectra are calculated without COBE-normalization (ar- 
bitrary units) and taking UHe ~ 4<r. (b) Matter power spectra 
for the same models as in (a). Normalization is arbitrary. 



On the other hand, increasing the baryon density has 
the well-known effect of increasing the height of the com- 
pression peaks (1st, 3rd,...) and decreasing the height 
of the decompression peaks (2nd, 4th,...). Of course, for 
very large cross sections, the effective viscosity of the 
CDM-baryon-photon plasma approaches that of a pure 
baryon-photon plasma and the CDM acts exactly like 
baryons. An effect similar to the increased viscosity due 
to CDM-baryon interactions can be seen in models with 
increased width of the last scattering surface due to non- 
standard recombination, where the effective viscosity of 
the baryon-photon plasma also increases around the time 
of recombination llTil . 
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The effect on the matter power spectrum seen in 
Fig. 0(b) can be understood in the same way. The in- 
creased viscosity in the interacting model leads to a much 
stronger damping on small scales than in the model with 
increased baryon density. 



III. COMPARISON WITH OBSERVATIONS 

In order to constrain models with CDM-baryon inter- 
actions we compare the CMB and matter power spectra 
of such models with recent observational data. 

CMB data set — Several data sets of high precision 
are now publicly available. In addition to the COBE |lS|| 
data for small I there are data from BOOMERANG jL9]f, 
MAXIMA |(]], DASI |fl| and several other experiments 
p2| . p3[ . Wang, Tegmark and Zaldarriaga [^2| have com- 
piled a combined data set from all these available data, 
including calibration errors. In the present work we use 
this compiled data set, which is both easy to use and 
includes all relevant present information. 

LSS data set — At present, by far the largest survey 
available is the 2dF of which about 147,000 galaxies 
have so far been analysed. Tegmark, Hamilton and Xu 
pBJ have calculated a power spectrum, P(fc), from this 
data, which we use in the present work. The 2dF data 
extends to very small scales where there are large effects 
of non-linearity. Since we only calculate linear power 
spectra, we use (in accordance with standard procedure) 
only data on scales larger than k = 0.2ft, Mpc~\ where 
effects of non-linearity should be minimal. 



The CMB fluctuations are usually described in terms 
of the power spectrum, which is again expressed in terms 
of Ci coefficients as 1(1 + l)Ci, where 



Ci = {\a lm \ 2 ). 



(28) 



The a; m coefficients are given in terms of the actual tem- 
perature fluctuations as 



(29) 



Given a set of experimental measurements, the likelihood 
function is 



£((-)') X cxp | ~x^[C{Q)- l ]x 



(30) 



where O = (f2, n&, H , n, r, . . . ) is a vector describing the 
given point in parameter space, a; is a vector containing 
all the data points, and C(O) is the data covariance ma- 
trix. This applies when the errors are Gaussian. If we 
also assume that the errors are uncorrelated, it can be 
reduced to the simple expression, £ cx e~ x ' , where 



N„ 



X 



E 

i=l 



(C/,obs — Ci, theory)! 



(31) 



is a x 2 -statistic and A max is the number of power spec- 
trum data points p9| . In the present letter we use equa- 
tion d3l| ) for calculating \ 2 . In the case where we also 
use LSS data % 2 is instead given by 



E 



,obs Cj, theory )^ 



+ 



E 



(P(fc)obs - P(fc)theory)j 



(32) 



r 



The procedure is then to calculate the likelihood func- 
tion over the space of cosmological parameters. The 2D 
likelihood function for (m c , a) is obtained by keeping 
(m c ,a) fixed and marginalising over all other parame- 
ters. 

As free parameters in the likelihood analysis we use 
f2 m , the matter density, f^, the baryon density, H , the 
Hubble parameter, n, the scalar spectral index, r, the op- 
tical depth to reionization, and Q, the overall normaliza- 
tion of the data. When large scale structure constraints 
are included we also use 6, the normalization of the mat- 
ter power spectrum, as a free parameter. This means that 
we treat Q and b as free and uncorrelated parameters. 
This is very conservative and eliminates any possible sys- 
tematics involved in determining the bias parameter. We 
constrain the analysis to flat (f2 m + Q\ = 1) models, and 
we assume that the tensor mode contribution is negligi- 



ble (T/S — 0). These assumptions are compatible with 
analyses of the present data [^2| , and relaxing them does 
not have a big effect on the final results. For maximiz- 
ing the likelihood function we use a simulated annealing 
method, as described in reference [p6| . 

Table I shows the different priors used. In the "CMB" 
prior the only important constraint is that 0.4 < h < 0.9 
(h = ff /(100 kms^Mpc' 1 )). For the CMB+H +BBN 
prior we use the constraint H$ — 72 ± 8 kms -1 Mpc -1 
from the HST Hubble key project |27j] (the constraint is 
added assuming a Gaussian distribution) and the con- 
straint fl b h 2 = 0.020 ± 0.002 from BBN @. Finally, 
in the iJo+BBN+LSS case, we add data from the 2dF 
survey pq) . 

Following Ref. we plot the constraints from var- 
ious experiments in the (m c , a) plane in Fig. 0; limits 
obtained with different priors are plotted with different 
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FIG. 2: Excluded regions in the (m c , a) plane for strongly- 
interacting dark matter, where m c is the mass of the dark 
matter particle, and a is the cross section for scattering be- 
tween dark matter and protons. Dashed and dotted lines are 
the CMB limits derived in this paper with various priors: dot- 
ted (red) is CMB only, short-dash (blue) is CMB+BBN+f/o, 
thick long-dash (black) is CMB+BBN+# +LSS. All of these 
limits assume a scattering cross-section between helium and 
dark matter of ane = 4<j. Thin dashed lines give the limits 
for CMB+BBN+# +LSS for spin-dependent scattering be- 
tween helium and dark matter (ane = 0, upper curve) and 
coherent scattering (a He = 16a, lower curve). In all of these 
cases, the region above the curve is excluded at the 95% con- 
fidence level. Solid lines give excluded regions from other 
papers. The region above the top solid line is excluded by 
the disk heating constraint from Starkman et al. Closed solid 
contours give excluded regions from the indicated space and 
balloon-based experiments discussed by Wandelt et al. The 
regions enclosed by these curves are excluded. The shaded 
(yellow) strip is the region suggested in Wandelt et al. for 
which self-interacting dark matter would have a comparable 
cross section for scattering off of both baryons and itself. 

line styles. The shaded region is the "Spergel-Steinhardt" 
region as given in equation ([!]); note that it is actually 
narrower than plotted in Fig. 3 of Ref. Q . This figure 
clearly shows that constraints from the CMB data alone 
for large CDM mass (m c > rrib ~ 1 GeV) are comparable 
to limits from galactic disk heating arguments jjj . On the 
other hand, for small CDM mass (m c < rrib ~ 1 GeV), 
the limit from the CMB alone is much stronger. This 
difference arises because the limits from disk heating are 
related to energy transfer, which decreases as (m c /m(,) 2 
for m c -C to;,, while the limits from the CMB and large 
scale structure are based on momentum transfer, which 
decreases only as m c /mb for m c <C nib- At the same time, 
the relative velocity also increases for low CDM mass (as 

m c ). This is the reason that, for m c <c mb, our bound 



goes as m c , but for the disk heating argument it goes 
as tti~ . Our bound is therefore much stronger at low 
mass, even if only CMB data is used. 

The addition of priors on and Hq does not increase 
the exclusion region by much, because there is very little 
degeneracy between (m c ,er) and these two parameters. 
Only when LSS data is added do we find a significant im- 
provement. The reason is that CDM-baryon interactions 
lead to a significant suppression of small scale power, 
as discussed in the previous section. Even for relatively 
small cross sections this effect leads to a disagreement 
with the 2dF data. With the inclusion of LSS data, our 
bound becomes much stronger than the disk heating ar- 
gument of Starkman et al. 0, even at large masses. A 
good analytical fit to our 95% C.L. limit is 

a 0.63a; 1 / 2 + 0.22x 3 / 2 

10- 24 cm 2 < l + a; 1 ^ (33) 

where x = m c /GeV. Equation ( |33| ) is accurate to within 
10%. For this case, we also consider the effect of coherent 
{o'He = 16it) and spin-dependent (er# e = 0) scattering 
between the dark matter and helium. The limits for these 
two cases are shown in Fig. 2 as thinner dashed curves 
above (for er# e = 0) and below (for one — 16a) the 
limit for a He — 4c. Altering a He hi this way changes 
the upper bounds by a factor of roughly 2 — 3, which is 
barely noticeable on the scale of this graph. 

IV. CONCLUSIONS 

CMB observations constrain any scattering interaction 
between dark matter and baryons. Our results indicate 
that the limit from current CMB observations is compa- 
rable to previous limits from disk heating for masses 
greater than 1 GeV, and the CMB limit is much stronger 
at smaller masses. If we also include large-scale struc- 
ture data, then our limit is more stringent that the disk- 
heating limit at all masses. An analytical fit of our com- 
bined CMB+LSS limit is given in Eq. (ff). 

Our CBM+LSS limits exclude the region discussed in 
Ref. in which self-interacting dark matter interacts 
with baryons with roughly the same cross section with 
which it interacts with itself (Fig. 2). This, by itself, 
does not exclude the self-interacting dark matter sce- 
nario, since the dark matter is not required to interact 
with baryons at all, but it does exclude models with the 
indicated dark matter-baryon scattering cross section. 

While space and balloon-based experiments can pro- 
vide tighter constraints for specific regions in the pa- 
rameter space of CDM mass and cross section, (Fig. 2 
and reference ||), the combination of CMB and large- 
scale structure appears to provide the best general upper 
limit on the CDM-baryon scattering cross section for ar- 
bitrary CDM masses (and it is much stronger than all 
other limits at low masses). Our constraints are some- 
what stronger for coherent scattering from helium nuclei 
and weaker for a spin-dependent interaction. Perhaps 
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more importantly, these limits will only get better with FG02-91ER40690. 
new CMB data. 
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TABLE I: The different priors used in the analysis. Notice that in the first two cases, the bias parameter b is not used because 
LSS data is not included in the fit. 

prior type fi m Q.bh 2 h n r Q b 

CMB Q 6 -l 0.008 - 0.040 0.4-1.0 0.66-1.34 0-1 free not used 

CMB + BBN + H Q b -1 0.020 ± 0.002 0.72 ± 0.08 0.66-1.34 0-1 free not used 

CMB + BBN + H + LSS Q b -1 0.020 ± 0.002 0.72 ± 0.08 0.66-1.34 0-1 free free 



